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ABSTRACT 



> 

■ It is believed that 77 Carinae is actually a massive binary system, with the 



wind- wind interaction responsible for the strong X-ray emission. Although the 
overall shape of the X-ray light curve can be explained by the high eccentricity 
of the binary orbit, other features like the asymmetry near periastron passage 
and the short quasi-periodic oscillations seen at those epochs, have not yet 
been accounted for. In this paper we explain these features assuming that 
, the rotation axis of 77 Carinae is not perpendicular to the orbital plane of 

^ . the binary system. As a consequence, the companion star will face ?/ Carinae 

on the orbital plane at different latitudes for different orbital phases and, 
since both the mass loss rate and the wind velocity are latitude dependent, 
they would produce the observed asymmetries in the X-ray flux. We were 
able to reproduce the main features of the X-ray light curve assuming that 
the rotation axis of 77 Carinae forms an angle of 29° ± 4° with the axis of 
the binary orbit. We also explained the short quasi-periodic oscillations by 
assuming nutation of the rotation axis, with amplitude of about 5° and period 
of about 22 days. The nutation parameters, as well as the precession of the 
apsis, with a period of about 274 years, are consistent with what is expected 
from the torques induced by the companion star. 

Key words: stars: individual: 77 Carinae - stars: binaries: general - stars: 
winds 
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1 INTRODUCTION 

The intensity and spectrum of the high energy X-ray flux, and its strict periodicity, are 
probably the strongest evidence of the binary nature of the rj Carinae system. The 2-10 
keV X-ray emission of rj Carinae is monitored by the Rossi X-Ray Timing Explorer RXTE 
since 1996, and the published results cover two cycles in the 5.52 year periodic light curve 
(Corcoran 2005). The duration of the shallow minima, as well as the general qualitative 
behavior of the light curve, were similar in the two cycles. The long lasting intervals of 
almost stationary intensity were modulated by low amplitude quasi-periodic flares, and the 
large flux increase that occurred before the minima was enhanced by strong short duration 
flares (Ishibashi et al. 1999). Although the X-ray light curve was successfully reproduced by 
analytical approximations involving wind-wind collisions (Ishibashi et al. 1999; Corcoran et 
al. 2001) and by numerical simulations (Pittard et al. 1998; Okazaki et al. 2008), which also 
reproduced the high resolution spectra obtained with Chandra (Pittard & Corcoran 2002), 
some features are still controversial, like the asymmetry near periastron passage, the short 
quasi-periodic oscillations seen at those epochs, and the difference in the phases of these 
oscillations between the two cycles (Okazaki et al. 2008; Parkin et al. 2009). 

Besides X-rays, other observational features can be related to wind-wind collision. Abra- 
ham & Falceta-Gongalves (2007) were able to reproduce the Hell A4686 line profiles and 
mean velocities detected close to the 2003.5 minimum by Steiner & Damineli (2004) and 
Martin, Davis & Koppelman (2006), assuming that they were formed in the cooling shocked 
material flowing along the winds contact surface. More important, to reproduce the line 
profiles reflected in the Homunculus polar cap (Stahl et al. 2005), they had to assume that 
the Homunculus axis is not perpendicular to the orbital plane. 

The rotation axis of rj Carinae probably coincides with the axis of the Homunculus; the 
shape of the nebula and the measured latitude dependent stellar wind velocity are strong 
indications that the rotational velocity is close to its critical value (Smith 2002; Dwarkadas 
& Owocki 2002). A consequence of the inclination of the rotation axis relative to the axis of 
the orbital plane is that the secondary star faces 77 Carinae at different latitudes as it moves 
along the orbit, and therefore, the latitude dependent velocity and mass loss rate of the 
primary's wind will affect the intensity of the X-rays produced in the wind-wind collision 
region. 
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Also, the large rotational velocity of 77 Carinae will affect its internal mass distribution, 
which will depart from spherical symmetry. The torque induced by the companion star 
will result in apsidal motions, as seen in other massive binary systems (e.g. see references 
in Claret & Gimenez 1993). Finally, the inclination of the rotation axis of 77 Carinae will 
produce nodding motions, which will further affect the strength of the wind-wind collision 
and the consequent X-ray intensity. 

In this paper we will take all these effects into account and calculate the X-ray light curve 
of 77 Carinae using the analytical approximation derived by Usov (1992) and the orbital 
parameters found by Abraham et al. (2005) and Abraham & Falceta-Gongalves (2007). We 
will show that for reasonable values of the precession and nutation periods and amplitudes, 
it is possible to reproduce the asymmetries in the light curve close to periastron passage and 
the amplitudes and phases of the short quasi-periodic oscillations for the two binary cycles 
observed by RXTE. 



2 THE X-RAY EMISSION MODEL 

We will use the model derived by Usov (1992) to calculate, at each point of the orbit, the 
X-ray luminosity originated in the shock heated gas at both sides of the contact surface, 
which depends on the mass loss rates (M p and M s ) and wind velocities (V p and V s ) of the 
primary and secondary stars, respectively, and on the distance D between them. Pittard & 
Stevens (2002) showed that for the 77 Carinae binary system, the major contribution to the 
X-ray flux comes from the interaction surface of the secondary wind, because of its higher 
temperature so that the expression derived by Usov (1992) and valid for adiabatic shocks 
becomes: 

^^^(^fw^'V->, (!) 



Aud 2 D V V s 

where d is the distance to 77 Carinae, taken as 2.3 kpc, and t(6 s ) the optical depth for X-ray 
absorption; M p and M s are expressed in units of 10~ 5 and 10~ 6 M yr _1 respectively, V p 
and V s in units of 10 3 km s -1 , and D in units of 10 13 cm; 9 S is the true anomaly, with 9 S = 
at periastron. We will not take into account any possible cooling of the very dense shocked 
gas very close to periastron passage (Parkin et al. 2009). 

We will assume that M s and V s have constant values, although Parkin et al. (2009) 
proposed a reduced secondary wind velocity near periastron to explain the observed change 
in the X-ray hardness ratio. On the other hand, we will assume that M p and V p depend on 
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orbital plane 





orbital plane 



Figure 1. Left: geometrical description of the intersection of the binary system orbital plane with the surface of r\ Carinae at 
latitude A, $ is the angle between the rotation axis and the perpendicular to the orbital plane; Right: definition of the different 
angles and coordinate systems involved in the precession of the line of apsis and nutation of the rotation axis of r\ Carinae 

the latitude A (9 S ) at which the orbital plane intercepts the side of rj Carinae that faces the 
secondary star, and can be expressed as (Dwarkadas & Owocki 2002): 

M(A) = M(90°)[l -fi 2 cos 2 A], (2) 
V(X) = V(90°) [1 -n 2 cos 2 A] 1/2 , (3) 

with f2 = uj/uj c ] to is the rotation velocity and uj c = (GM P / Rp) 1 ^ 2 its critical value; G is 
the gravitational constant, M p and R p are the mass and radius of rj Carinae, respectively. 
These expressions are valid when Q is close to unity; they were already used to reproduce 
the observed wind velocity as a function of latitude in rj Carinae, as well as the shape of the 
Homunculus nebula (Smith 2002; Dwarkadas & Owocki 2002). 
By replacing eq. (2) and (3) in (1) we obtain: 



F x {0. 
with 



G(t) 



M p {\) 



D LM P (90°) 



1/2 



VJ9Q ) 



1/2 
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M p (90 o )^ p (90 o ) 
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(4) 



(5) 



where we have defined t(9 s ) = t p (9 s ) + r (t); r p {6 s ) represents the absorption produced by 
the wind of 77 Carinae intercepting the line of sight, and To(t) is constant or a slowly varying, 
phase independent function of time, representing all other sources of absorption. 



2.1 The effects of precession and nutation on X(9 S ) 

r) Carinae must be highly distorted, both by rotation and by the presence of the companion 
star in a highly eccentric orbit. Although the total angular momentum in a detached binary 
system is conserved (except for a small amount lost by the stellar winds), energy will be 
dissipated by the tidal forces, until a minimum energy equilibrium configuration is reached, 
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in which the orbit is circular, and the stellar spins perpendicular to the orbital plane and 
equal to the orbital velocity. However, although the orbital parameters are continuously 
changing, the orbit can be consider instantaneously Keplerian, and precession and nutation 
rates can be calculated as averages over the orbital period (Eggleton, Kiseleva & Hut 1998). 

In this section we will derive the effects of precession and nutation on the latitude \{6 S ). 
We will assume that n Carinae rotates with angular velocity u around an axis that forms an 
angle $ with the perpendicular to the orbital plane [z axis in Figure 1), and its projection 
on the orbital plane forms an angle with the line of apsis (x axis), so that: 

sinA(# s ,t) =sine(0 s ,t)sin$(t), (6) 

where 

tan0(M) = — , (7) 



sin $(t) = AL^ ». ( 8 ) 

UJ 

As a consequence of nutation, the rotation angular velocity vector ul will describe a cone of 
amplitude Aip and period P n around its non-perturbed direction, which forms an angle ip 
with the polar axis z (Figure 1). 

In a coordinate system (x", y", z"), in which z" is directed along the unperturbed rotation 
axis, the components of ul will be (a; sin Aipsmuj n (t), a; sin Atp cosuj n (t), cjcosAy?), with 
oj n (i) = (2-n/P n )At + 9 0n ; 9 0n is a constant phase and At = t — t , with t = 2,450,795 
being the JD of the beginning of the 1997.9 minimum. To obtain the components of ul in the 
(x,y,z) coordinate system, in which z coincides with the orbital axis, the (x",y",z") axis 
must be rotated around y" by an angle </?, so that the new z' axis coincides with z, and then 
around z' by and angle 6' s = 9 S + 2n/ P p + 9 0p , to take into account the orbital motion and the 
precession of the line of apsis; 9q p is a constant phase. These rotations can be represented 
by the matrix M: 



M = 



' cos9' s sm9' s 
-sin^ cos9' s 



^ / COS (f 



\ 







I J 



— sin ip 



sin<^ 
cosy? 



(9) 



V 

so that (x, y, z) T = M(x", y", z") T 







/ 
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Figure 2. Observed 2-10 keV X-ray flux, from Corcoran(2005), shown as open circles, and model, shown as a continuous line, 
obtained with the parameters listed in Table 1, but with Aip = 0° (no nutation) in the two upper graphs and Aip = 4?5 for 
the lower ones. 
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2.2 Opacity and the Hi column density 

As mentioned before, the optical depth for X-ray absorption was divided into two parts: 
r(9 s ) = T (t) + r p (9 s ), where r p (9 s ) = o- ph N H (9 s ) represents the photoelectric absorption 
produced by the unshocked wind of rj Carinae intercepting the line of sight to the vertex of 
the X-ray emitting cone, with column density Nh, and 7o(£) that represents all other sources 
of absorption, excluding the stellar wind; cr ph is the cross section of photoelectric absorption, 
multiplied by the heavy element abundance relative to H. N H (9 S ) is calculated from: 

Nh{9 , )= _l_ &\r * (10 ) 

4-7rfmi H \V P ) Js Bh s 2 + rfi 

where fi is the molecular weight and m H the mass of the hydrogen atom, M p and V p are 
the mean values of the mass loss rate and wind velocity of the primary star; s is measured 
along the line of sight to the apex of the X-ray source; r = 6 sin \1> and s sh = b cos *, where 
b is the distance from rj Carinae to the shock, measured in the orbital plane; \& is calculated 
from: 

sin\l/ = sin(# s — 9 ) sini, (11) 

where i is the inclination of the orbit and 9 is the true anomaly at conjunction. While 
the inclination is one of our model parameters, the value of the input parameter 9 is still 
controversial (Pittard et al. 1998; Corcoran et al. 2001; Falceta-Gongalves, Jatenco-Pereira & 
Abraham 2005; Kashi & Soker 2007; Hamaguchi et al. 2007; Abraham & Falceta-Gongalves 
2007; Okazaki et al. 2008; Falceta-Gongalves & Abraham 2009; Parkin et al. 2009). As we 
will see later, it affects mostly the opacity near periastron, where the model anyway fails to 
reproduce the duration of the shallow minima. 

By solving the integral of equation (10), we can write: 

tJ9 s ) = ; ^ T (— — arctan — ^-—] , (12) 
pV ; 6sin^ V2 tan*/' v ; 

where C T = a v ^M p j^\imHV v - 

The contribution to the opacity of the unshocked wind of the secondary secondary star is 
much smaller than that of the wind of i] Carinae, because of the much smaller value M s /V s . 
However, depending on the position of the secondary star on the orbit near periastron, 
the absorption due to the shocked gas intercepting the line of sight could be large (Falceta- 
Gongalves, Jatenco-Pereira & Abraham 2005; Parkin et al. 2009) and can affect the duration 
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Table 1. Input and model parameters 



Input Model 



e = 0.95 P p = 274 ± 15 years 

9 = -45° 9 0p = 5?3 ± 1° 

P = 2024 days P n = 22.45 ± 0.04 days 

A = 15 A.U. 9 0n = 90° ± 10° 

t = 2, 450, 795 JD H = 0.975 ± 0.001 

ip = 29° ± 4° 

A(p = 4?5 ± 0?5 

C x = (7 ± 1) x 10~ 5 

C T = 7.7 ±0.5 A.U. 

Go = (2200 ± 150) X 10~ n erg cm"^- 1 
i = 60° 



Table 2. Inclination of the binary system orbit i and of rj Carinae rotation axis i* 



i i*(r) = 0.2) i*(ri = 0.1) 
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of the minima in the light curve; this issue will not be addressed here since it requires 
numerical simulations. 



2.3 The orbital parameters 

We will use the orbital parameters derived by Abraham et al. (2005) from the observed 
7-mm light curve of 77 Carinae during the 2003.5 minimum and listed in the left column 
of Table 1. They were successfully used to reproduce the HenA4686 line profiles and their 
mean velocities, although to reproduce the mean velocities and profiles of the emission lines 
reflected in the polar cap of the Homunculus, it was necessary to assume that the Homunculus 
axis is not perpendicular to the orbital plane. Considering that the angle between this axis 
and the line of sight is 45° (Davidson et al. 2001; Smith 2006), for each inclination % of the 
orbit relative to the observer, an orientation for the Homunculus axis i* was found, for which 
the reflected line profiles and velocities could be reproduced (Abraham & Falceta-Gongalves 
2007). In Table 2 we present the values of these angles for two values of rj = M S V S /M P V P . 
In the next section we will use the value of i* — tp + Atp, obtained from the model that 
reproduced the observed X-ray light curve, to constrain the value of the orbital inclination 
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Figure 3. Details of the X-ray emission model without precession and nutation, close to the epochs of the two shallow minima, 
for three different values for the precession phase 8op- —4°, 5° and 13°. All other model parameters are shown in Table 1. 




-180 -135 -90 -45 45 90 135 180 
Orbital Phase (degrees) 

Figure 4. Variations with orbital phase of the wind velocity, mass loss rate, wind density at the contact surface, and latitude 
at the intersection of r\ Carinae with the orbital plane facing the companion star. The first two quantities, given by the left axis, 
are relative to their values at A = 90°, the latitude values, in degrees, arc given in the right axis. The wind density, relative to 
its value at A = 90°, in units of 1.5 X 10 11 cm -1 is also displayed at the right axis. 



3 RESULTS 

We used equations (1) to (12) to model the X-ray light curve of r\ Carinae, with the orbital 
parameters listed in the first column of Table 1. No formal fitting was attempted; instead, 
the model parameters were changed until they reproduced the general shape of the light 
curve (except the shallow minima), and the amplitude and period of the oscillations that 
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occurred just before them. The parameters that fulfilled these criteria are listed in the last 
column of Table 1, and the model superposed to the observed light curve is presented in 
Figure 2; the two upper graphs represent the model without nodding motions (Ay? = 0), and 
the last two graphs include nodding. The multiplicative term G(t) in equation (5) was fitted 
by a function G(t) = G [l + C x (t — t )], with G and C x having constant values. The quoted 
errors were obtained by changing the values of each parameter while keeping the others 
constant until the model was no longer acceptable. As mentioned before, no formal fitting 
was attempted, but the combination of model parameters could not be changed arbitrarily, 
since each of them affect some particular feature of the X-ray light curve, as discussed bellow. 

The parameters Q and 9 0p are responsible for the asymmetry in the light curve at both 
sides of the shallow minima, as can be seen in the first two graphs in Figure 2. To illustrate 
this dependence, we plotted in Figure 3 the model light curves close to the minima for 
Q = 0.975 and several values of 6q p , when no precession or nutation are present. We found 
that 6q p = 5?3 reproduces well the asymmetry in the mean value of the light curve close to 
the minimum of 1997.9, but 9o p ~ 13° was necessary to reproduce the minimum of 2003.5. 
The difference in angles was attributed to the motion of the apsis, resulting in the precession 
period listed in Table 1. A large value of was necessary to get the observed discontinuity; 
it indicates that 77 Carinae is rotating at almost its critical velocity, as expected from the 
episodes of large mass loss. Of course, if other causes were responsible for the asymmetry in 
the X-ray light curve, like changes in opacity or stellar wind parameters, the value of these 
parameters should be revised. 

The nutation parameters P n , 9 0n and Aip are responsible for the amplitude and period of 
the large oscillations that occur before the minima, as can be seen in the two lower graphs in 
Figure 2; it is important to notice that, near periastron, the oscillations remained in phase 
in the two cycles for the same initial phase 8q h . However, the model does not reproduce the 
amplitude and the period of the oscillations far from periastron, which could be expected 
from the variation of the torque of the secondary star acting on i] Carinae. 

The opacity parameter C T determines the shape of the light curve before the minima, 
and has no influence at other phases. The linear dependence of G(t) with time implies an 
increase of about 30% in the X-ray intensity between the two cycles, which coincides with 
an increase in the optical flux during the same time interval (Martin et al. 2004) and could 
maybe attributed to an overall decrease in opacity. 

No assumptions were made on the values of the mass loss rates and wind velocities of the 
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binary stars. Instead, the value obtained for the model parameter Gq and the comparison 
between the model hydrogen column density Nh and that inferred from the XMM X-ray 
spectra observed far from periastron passage (Hamaguchi et al. 2007) allowed us to put some 
constrains on their magnitudes. From eq. (5) and assuming r = we can write: 

M 2 G 



s 

V, 1.83 xlO- 10 ' 



^ 1/2 (90°) (13) 

l.OO X 1U iv 

where 

„ (90 °) = . ^ = mm „ (A) (14) 

M p (90°)Vp(90°) M p (90°)\/ P (90°) 

From our model r?(90°) = 0.15^(30°), resulting that for n{30°) = 0.2 and V s = 3000 km 
s -1 , M s = 8 x 10~ 6 M Q y" 1 , well within the values of the mass loss rate of the secondary 
star used in the literature. 

The hydrogen column density inferred from the X-ray spectra in January 2003 was 
Nh = 9 x 10 22 cm~ 2 (Hamaguchi et al. 2007). Using this value in equation (10) we obtain 
Mp/V p = 4.4 x 10 14 g cm" 1 , where we have used: 



is. 



" s [ /JT - arc t ail -J— ) = 0.062 AT (15) 



s sh s 2 + rl 6sin^\2 ' tan*. 

Assuming a wind velocity of 500 km s" 1 for rj Carinae, we find M p = 3.3 x 10" 4 M y _1 , 
also consistent with the values found in the literature. 

Finally, the model opacity for January 2003 was r p = 0.69, which together with the 
observed hydrogen column density gives a value for a p h = t p /Nh = 4.9 cm 2 g" 1 , consistent 
with the opacity to 3 keV photons of a 10 4 -10 6 K gas (Parkin & Pittard 2008). 

From the inclination of the rotation axis of i] Carinae we estimated that, relative to the 
observer, the orbit has an inclination i ~ 45° — 60° , for rj varying between 0.2 and 0.1, as 
can be seen in Table 2. The variation of the wind velocity, mass loss rate and wind density 
along the orbital period, as well as the variation in A are shown in Figure 4 for the first 
orbital period. 

Although the eccentricity and period of the binary orbit used in the model agree with 
those used by other authors (Pittard et al. 1998; Ishibashi et al. 1999; Corcoran et al. 2001; 
Pittard & Corcoran 2002; Okazaki et al. 2008), the value of 9 is still controversial. The model 
presented in this paper was calculated for 6q = —45°; changing its value to #o — +45° will not 
have any effect on the X-ray emission but affect the absorption, mostly close to periastrom 
passage, where the model does not anyway reproduce the X-ray light curve; however, using 
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Figure 5. Asymmetry parameter k of the internal mass distribution of r) Carinac , produced by rotation and gravitational 
torque induced by the companion star, versus the ratio between the primary star radius and the orbital semi-major axis, for 
three values of the primary to secondary mass ratio M s /M p . Horizontal lines are the expected values of k for polytropes of 
index n = 0, 1, 2 and 3. 

#o = +45° will not reproduce the shape and central velocity of the HenA4686 and Paschen 
lines (Abraham & Falceta-Gongalves 2007; Falceta-Gongalves & Abraham 2009). 



4 DISCUSSION 

As mentioned before, the model does not adjust the amplitude and period of the oscillations 
in the X-ray light curve far from periastrom passage. From Figure 2 we can also see that 
it does not reproduce the duration of the shallow minima. In fact, neither the analytical 
models nor the numerical simulations developed up to the present time were able to account 
for the extended minima as a result of X-ray photoelectric absorption by the dense wind of 
i] Carinae intercepting the line of sight (Pittard et al. 1998; Ishibashi et al. 1999; Corcoran 
et al. 2001; Pittard & Corcoran 2002; Hamaguchi et al. 2007). Possible explanation are the 
increase in the H column density to ~ iq 24-26 cm -2 due to additional material provided 
either by a slowly expanding shell of shocked material formed during periastron passage 
(Falceta-Gongalves, Jatenco-Pereira & Abraham 2005), or by the primary wind itself, which 
"engulfs" the secondary star (Okazaki et al. 2008), or simply by the suppression of the 
secondary wind due to accretion of matter from the close primary star (Soker 2005; Akashi, 
Socker & Behar 2006) . 
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Figure 6. Angle A(p between the total and orbital angular momenta as a function of R p /A for several values of the mass ratio 
M a /Mp. The horizontal line shows the value of A(p obtained from the model. 
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Figure 7. Time-scales for tidal friction txF, alignment t a n g , circularization t c j rc , and synchronization t sp i n for (M s /M v = 0.4). 
The vertical line separates the regions where the orbit is unstable (left) and stable (right). 

4.1 The orbital stability 

The precession and nutation parameters are related to the torques of the secondary acting on 
the fast rotating non-spherical primary star and can provide constrains on the stellar masses, 
internal structure and orbital stability. Hut (1982) derived expressions for the evolution of 
the orbital and spin parameters in highly eccentric binary system, assuming that the tidal 
bulges lag by a constant angle from the line that joins the stars. Eggleton et al. (1998) 
obtained an expression for this angle, considering that the dissipative force is proportional 
to the rate of change of the quadrupole tensor of the stars, as seen by an observer that 
rotates with them. 

The precession rate of the apsis, is independent of dissipation and, neglecting the torques 
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Figure 8. Relation between the mass ratio M s /M p and R p /A, for the combination of model parameters (<p, Aip) = (34°, 5?0), 
(29°, 4? 5), and (24°, 4? 0). The shadowed region represents the parameter space in which the orbit is unstable. 



of 7] Carinae on the secondary star, can be expressed as: 
2 (3 cos 2 up - 1) fR p 



with 



/(e) = (l-e 2 )- 5 (l + ^ 2 + ^e 4 ) (17) 
9(e) = (1-e 2 )- 2 , (18) 

where k is the constant part of the quadrupole moment, which depends on the internal mass 
distribution of the primary star, A is the semi-major axis of the orbit, and R p is the radius 
of the primary star. 

In massive binary systems, the measured precession period P p , together with the orbital 
parameters, stellar masses and radii, are used to calculate k and improve the stellar structure 
models (e.g. Claret & Gimenez 1993). Since for 77 Carinae the orbital and stellar parameters 
are unknown, we will take M s /M p and R p /A as free parameters. 

In Figure 5 we show the values of k that satisfy equation (13) for P/P p = 0.02 (obtained 
from our model) as a function of Rp/A, for several values of M s /M p . The maximum value 
allowed for R p /A corresponds to the separation between the stars at periastron: R v j A = 
(1 — e) = 0.05. The horizontal lines represent the values of k for rotating polytropes with 
indices n — 0, 1, 2 and 3 (Chandrasekhar 1933); n = corresponds to a star with constant 
density while for n — 1 the radius is independent of the central density. 

Another constrain for the ratios M s /M p and R p /A can be obtained from the amplitude 
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of nutation Atp, which represents the angle between the orbital and total angular momenta, 
and can be expressed as: 

cot Atp = cot tp + - — - , (19) 

Ioj sin tp 

where h is the orbital angular momentum, I = M p r 2 g R 2 p the momentum of inertia of the 
primary star and r 2 a parameter that depends on its internal structure. We can also write 



h/Iuo in terms of M s /M p and R p /A: 
h _ (1 + MjMp) 1 / 2 (RpX 1 ' 2 



(20) 



Iuj M s /M p \AJ r 2 g £l(l - e 2 ) 1 / 2 ' 

In Figure 6 we present the relation between Atp and R p /A for tp = 29° (obtained from 
our model), and M s /M p = 0.1, 0.4, and 1.0. In the calculation, we used an interpolated 
relation between r 2 and k for rotating polytropes obtained from Motz (1952): 

log r 2 = 0.453 log k - 0.307 (21) 

The value of Atp = 4.5 derived by fitting our model to the X-ray data is also shown in the 
figure. We can see that M s /M p ~ 0.4 represents the minimum value of the mass ratio for 
which a solution can be found. However, this result depends of the actual value of r 2 , and 
should be consider only in the context of a consistency test for the parameters derived from 
the X-ray light curve. 

The last parameter derived from the observations is the nutation period, which according 
to Eggleton et al. (1998) depends on both the conservative torques and those produced by 
tidal dissipation and can be written as: 



P k M s (^ M s \- X l 2 (R p \ z l 2 cosy? 



P n 2r 2 M p { 1 + M p ) [a) (1-e 2 ) 3 / 2 (22) 
3 M s f M s y l ( A \ 2 e 2 (l + l/6e 2 ) P 



8r 2 g M p \ M P J \R P J (1-e 2 ) 9 / 2 2vrt TF ' 
where £tf is the tidal friction time scale. 

When only non-dissipative torques, represented by the first tem in eq. (19), are con- 
sidered, the nutation period turns out to be several orders of magnitude larger than that 
obtained from our model (P n = P/92) for any combination of k, (Rp/A) and (M s /M p ) 
given by equation (13); therefore, the observed nutation must be produced by the dissi- 
pative torques, i.e. the second term in eq. (19), and its value can be used to estimate the 
dissipative time-scale trF- 
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3 M s 



8r 2 M p 



M, 



A\ 2 e 2 (l + l/6e 2 )P„ 



1?, 



(23) 



(l_ e 2)9/2 2 7t' 

In Figure 7 we display as a function of R p /A for M s /M p = 0.4. 

When the tidal dissipation time-scale is known, it is possible to determine the time-scales 
for spin alignment and synchronization, and orbit circularization, defined as: 



y_ 



dy 



dt 



(24) 



where t y = (t alig ,t spin ,t circ ) for y = (<p,u,e), respectively. 

For the case of a highly eccentric orbit, they can be obtained from Hut (1982): 
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In Figure 7 we plotted t a u g ,t sp i n , and t C i rc as a function of R p /A for M s /M p = 0.4. 
We can see that although Ttf is of the order of 10 5 years, the timescales for alignment, 
synchronization and circularization are larger than 10 9 years, which means that the binary 
system did not reached yet its equilibrium configuration, considering an evolution timescale 
of 10 6 — 10 7 years for the massive stars. 

We should notice in eq. (27) that H 3 can be negative, implying that the eccentricity can 
increase with time, which corresponds to an unstable orbit. This occurs when the spin of 
the primary star is larger than the orbital angular velocity at periastron, or: 

5> < l ~ e2 1 (37) 

A (1 + e) 3 /2 (l + M S /M P )V3 ^ 0> > 

For M s /M p = 0.4 this occurs when R p /A = 0.035, this value is shown as a vertical line 

in fig. (7). The general relation given by eq. (34) is shown in Figure 8; the shadowed area 

represents the parameter space for which the orbit is unstable. In the same Figure we show 

the lines that satisfy eq. (16), for the parameters found in our model (</? = 29°, Acp = 4?5), 

and for their extremes given by the uncertainties (tp = 34°, Atp = 5?0) and (tp = 24°, A(p = 

4?0). The values of M s /M p and Rp/A for the binary system should lie between these two 

extreme lines. 



5 CONCLUSIONS 

We were able to reproduce the general features of the 2-10 keV X-ray light curve of 77 Carinae 
obtained by RXTE (Corcoran 2005) during two cycles, including the amplitudes and phases 
of the short period oscillations that occurred prior to the shallow minima, assuming that 
the star rotates around an axis that is not perpendicular to the orbital plane, so that the 
secondary star faces 77 Carinae at different latitudes as it moves along the orbit, using the 
fact that both the mass loss rate and the terminal wind velocity of 77 Carinae are latitude 
dependent. According to the model, the star should be rotating with a fraction Q = 0.975 
of its critical velocity around an axis that forms an angle of about 29° with the axis of the 
orbital plane, nutates with an amplitude of about 5° and a period of 22.5 days. We also 
found that the line of apsis precesses with a period of about 274 years. According to the 
results of Abraham & Falceta-Gongalves (2007), the inclination obtained for the rotation 
axis of 77 Carinae implies that the inclination of the binary orbit relative to the observer 
must be % ~ 45° — 60°, depending on the ratio of the wind momenta (77 ~ 0.2 — 0.1). 
To fit the overall amplitude of the X-ray light curve, we calculated the phase and latitude 
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dependent opacity due to the wind of rj Carinae intercepting the line of sight, and introduced 
a phase independent absorption, which decreased linearly with time and explained the overall 
increase in the X-ray flux between the two cycles. We used the precession period and nutation 
amplitude and rate to constrain the mass ratio of the binary system and the radius of the 
primary star relative to the semi-major orbital axis. We found that the orbit is stable if the 
radius of rj Carinae is larger than 0.035 times the orbit major axis, and the mass of the 
companion star at least half the mass of 77 Carinae. Finally we found that for stable orbits, 
the time scale for orbit circularization, spin alignment and synchronization is much larger 
than the lifetime of the stars. 
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